On the evolution of cosmic-superstring networks 
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We model the behaviour of a network of interacting (p, q) strings from IIB string theory by 
considering a field theory containing multiple species of string, allowing us to study the effect of 
non-intercommuting events due to two different species crossing each other. This then has the 
potential for a string dominated Universe with the network becoming so tangled that it freezes. We 
give numerical evidence, explained by a one-scale model, that such freezing does not take place, 
with the network reaching a scaling limit where its density relative to the background increases with 
N, the number of string types. 
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I. INTRODUCTION 



Despite the remarkable success of inflation 0, Q 
in describing the observed patterns of the microwave 
anisotropy .3] there remains an active interest in cos- 
mic strings as a serious player in the early universe 0. 
Whilst it is now clear that cosmic strings cannot play a 
dominant role in structure formation their appearance 
in many scenarios, whether they be motivated by M- 
theory or supergravity leads us to consider them as a 
sub-dominant partner 0, IE with tensions of the order 
Gfi < 1CP 6 . Depending on the model and the allowed de- 
cay processes for the string network, the bounds on Gp, 
can be significantly strengthened. These could include 
the production of dilatons 0, of gravitational waves 
or even the possibility of loops of string being trapped 
in the compact dimensions leading to a monopole type 
problem for the strings [lfj. 

For many years the main connection between cosmic 
strings and superstrings was that both were much studied 
but never observed. Indeed, the high tension of funda- 
mental strings, placed at around the Planck scale, ruled 
then out as cosmic string candidates. However, over the 
past decade or so, more interesting solutions have begun 
to emerge in string theory, involving warped compactifi- 
cations or large extra dimensions. One of the by-products 
has been a realization that these warp factors play a cru- 
cial role in determining energy scales of compactincations 
In string theory these warp factors are sourced by 
brane-fluxes and much has been made of their use in 
model building[ll[Ii. 

Recently the possibility that the same basic mechanism 
could provide both the inflaton and lead to the formation 
of cosmic strings has received renewed attention in the 
context of braneworlds E E 0, E 0. Some 
consistency conditions for superstrings to act as cosmic 
strings were given in 0, in particular the KKLMMT0 
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picture of IIB string theory has a warp factor which al- 
lows the Dirichlet-fundamental strings (p, q) bound state 
to have reduced tensions and behave as cosmic strings. 
Such cosmic strings are different to the standard Abelian 
strings as they offer more ways of interacting; two differ- 
ent (p, q) strings could pass though each other, or they 
could reconnect in one of the two ways depicted in Fig. 
^20] . A network formed from such strings will therefore 
have quite different properties. Networks formed from 
Abelian strings consist of loops and long strings, (p, q) 
networks also contain loops and long string, but there 
are also links which start and end at a three-point vertex 
Fig. H The presence of these links in non- Abelian string 
networks has lead authors to speculate the possible exis- 
tence of a frozen network which comes to dominate the 
matter content of the Universe 0, |2^, A number 

of authors have begun to reconsider the evolution of the 
more complicated dynamics associated with cosmic su- 
perstrings, and arg ued in general that scaling solutions 
are achieved pil l25l l2^| . For example, recently there has 
been work suggesting that such networks of (p, q) strings 
reach a scaling behaviour, where the density of strings 
scales with that of the background energy density of the 
Universe • Here we re-consider the field theory model 
of j^, which can be thought of as containing N basic 
strings and ^N(N — 3) bound states to give a total of 
^N(N — 1) different species of string. We shall see from 
the numerical simulations that as N increases, so does 
the scaling density. This is to be expected because the 
network loses string via the formation of loops, and loops 
can only form from the reconnection event caused by two 
strings of the same type intersecting. Because the chance 
of a string meeting the same type decreases for large N 
then the probability of loop formation at reconnection 
goes down. Although we do not claim this model accu- 
rately describes a network of (p, q) strings it does have the 
important property that strings reconnect as they pass 
through each other, so creating the links which are the 
main difference to the standard cosmic string scenario. 
An important difference between our model and that of 
the (p, q) relates to the allowed string tension. While our 
field theory gives all strings the same tension, the tension 
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of a (p, q) string increases for larger p, q having the effect 
that low p, q will be preferred dynamically. From the per- 
spective of the field theory model this would correspond 
to having N small. 

Before discussing the dynamics associated with a net- 
work of non-abelian strings, we first remind the reader 
about some of the basic properties associated with or- 
dinary abelian strings. For a comprehensive review see 

The initial distribution of a network of strings is com- 
prised of a combination of infinite strings that stretch 
across the universe and a scale invariant distribution of 
loops of string. Roughly 80% of the string is initially 
infinite in length and 20% is in loops. There exist a se- 
ries of characteristic length scales on the network, but 
the one we are interested here is £, which corresponds 
to the typical separation between strings, and is given in 
terms of the volume V and total length L of strings by 
£ 2 = V/L. Other possible scales correspond to the typ- 
ical curvature scale on the strings, and the small scale 
structure on the string network. As the universe ex- 
pands, the network evolves and the characteristic length 
scales change with time as well. There are a number 
of physical processes that occur allowing the network to 
lose energy. The long strings can cross themselves form- 
ing loops which are chopped off the long strings, strings 
oscillate and because they are so massive they radiate en- 
ergy through gravitational radiation. For our case, where 
we are considering a field theory model, the prominent 
decay mode will be through particle production. Mean- 
while the long strings are stretched by the expansion, 
which increases the energy stored in the network. Hence 
there are competing effects which act against each other. 
Determining precisely how such a network evolves is a 
challenge that has proved difficult to solve exactly. The 
equations of motion are non-linear and the requirement 
of keeping track of the whole network in case new loops 
are formed from string intersections, requires consider- 
able computing power. The net result is that the correla- 
tion length £ scales with the Hubble radius, leading to the 
ratio of the string energy density and the background ra- 
diation/matter energy density becoming a constant, im- 
plying that the strings provide a fixed contribution to the 
total energy density in the universe. This scaling solution 
formed the backbone of investigations into string dynam- 
ics. The initial loop distribution is modified to one where 
the number density of loops of size I falls off as a power 
law in both the radiation and matter dominated eras, 
whilst maintaining the scale invariant distribution. The 
case for non-abelian strings is not as well understood, not 
least because much less effort has gone into following the 
evolution of these more complicated objects. Aryal et al 
[29| investigated the formation of Z$ strings. Based on 
numerical simulations they argued that unlike the case 
for Abelian strings, the initial loop distribution was not 
a scale invariant power law distribution but rather it was 
better approximated by an exponential distribution with 
the number density falling off exponentially with loop 





FIG. 1: Example of formation of string segments due to non- 
intercommuting . 



length. The majority of the string length is initially to 
be found in one network taking up 93% of the total string 
length, with the rest in loops and much smaller networks. 
It raised the important question, what happens to such 
a non-abelian network, and this was partly addressed in 
Refs. |22, H3, 0n|? In this paper we will look more 
closely at the evolution of non-abelan strings, and show 
that they always appear to enter a scaling regime with a 
corresponding loop distribution that, rather than having 
a scale invariant power-law distribution, instead appears 
to be exponential. We will follow and modify the ap- 
proach developed by Vachaspati and Vilenkin in Ref. [3(j 
where the authors investigated the evolution of a network 
of strings connected to monopoles. 



II. ONE SCALE MODEL 

In this section we shall describe an analytic approach 
to understanding the scaling regime of a network of non- 
Abelian string. In what follows we will assume that the 
network evolution is described by a single length-scale, 
£, and shall see that the solution to this model is that 
a network evolves in a self-similar fashion with £ ~ t, 
just as for Abelian networks. We start by writing down 
the evolution equation for a network containing a single 
string species evolving in an expanding background with 
Hubble parameter H — a/ a, where a is the scale factor 
and a 



dt ' 



Taking the pressure of the string network 
as p — 7/3, where p is the energy density of the strings, 
and 7 is the equation of state parameter, then the fluid 
equation is 



(2.1) 



The first term on the right hand side includes the work 
done by the system per unit volume per unit time, and 
the second term describes t he p rocess of energy loss from 
strings by particle emission [^3] • Taking p to be the string 
tension we define the network length-scale £ through the 
string density, 
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For a scale factor a oc t v , then 
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where k = |(1 + 7)^, A 



% . This can be solved to give 



III. NUMERICAL MODEL 



m = 



1 



-t + Kt K . 



(2.5) 



Causal evolution of the network requires £ to grow no 
faster than £ so we must have k < 1, leading to £ oc £ 
asymptotically, typical of self-similar evolution. 

Now we include loop production and reconnection. 
Following 30], normally we would expect that most of 
the loops chopped off the network are likely to be re- 
captured, both events occurring on a time-scale of order 
£(£). However, including expansion, the two events are 
not quite equal. In a time interval £, the density of strings 
is reduced by a fraction £/£, leading to a similar reduc- 
tion in the probability of loop absorption. In other words 
there is an imbalance with the infinite network losing en- 
ergy to loops at the rate pi 00 ps ~ — (£A)(m£/£ 4 ) ° 1 ' 



Pioops ~ -cip/t, 



(2.6) 



where we expect the constant q ~ 1. Things change 
even more in our particular case. We shall see later that 
there is not one type of string but ^N(N — 1), which 
can be thought of as bound states of N types of under- 
lying string. Now, imagine two of these strings Ni and 
N2 which cross one another. Only when N± and N% are 
the same will they form loops. Hence the probability of 
forming loops is decreased in our case by an additional 
factor of 2/(§iV(iV - 1) + 1). We see this by noting that 
for n(— t}N(N — 1)) string species there are \n(n + 1) 
distinct pairings of string, of which n will pair up the 
same type. So the chance of a random pairing contain- 
ing two of the same species is 1 . = 1 , r , Ar 2 , , In 

other words we have 

2ci 



^JV(JV-l)+l 



^N(N - 1) + 1 £ ■ 
In the rate equation l|2.1|l we now have 



p = -3H(l + 7 )p-c pP /^ 
= -2np/t - c p p/£ 
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(2.8) 
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Note the form of the loop term is simply to re-normalise 
the work-done contribution. In other words the solution 
for the characteristic length scale is as we quoted earlier 
but with a new k: 



X 



1 — Kpj 



£ + K N t^ 



(2.10) 



-. Now we have obtained the 



where k n = k + t^ 1)+1 . 
characteristic length of the network as a function of N 
we can compare this to our numerical solutions, which 
we now describe. 



In order to model a network of strings which do not 
intercommute, but rather form links between the two 
string segments as they cross (see Fig. ^| we followed 
the techniques described in 22] for setting up a linear 
sigma model. This is an approximation for describing a 
scalar field in the long wavelength limit where the field 
is taken to live on its vacuum manifold. In order to have 
different types of string this manifold must have more 
than one non-contractible loop, with each such loop cor- 
responding to a string solution. We take the vacuum as 
given in [22| , and depicted in Fig. [21 which takes the form 
of N spokes whose ends are identified, noting that N = 2 
corresponds to the usual case of a single species of string 
[33|] . In general there are ^N(N — 1) species correspond- 
ing to the number of incontractible loops in the vacuum. 
For example, referring to Fig. [21 we could form strings 
from the loops AOBA, AOCA, BOCB,... Therefore we 
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FIG. 2: Vacuum manifold for the sigma model, the ends of 
the spokes are identified. 

evolve a scalar field which takes values in [0, ir) x Zjv and 
has a canonical kinetic term. The initial conditions are 
chosen such that the field takes random values at each 
lattice site and we use a 300 3 lattice. The first set of 
simulations were performed in Minkowski space and the 
results for the network length-scale arc given in Fig. [31 
showing how £ evolves in time for various N. Note that 
the network approaches a scaling solution in each case, 
with £ oc £. Also, in the simulations with a larger num- 
ber of species the network grows more slowly in time, as 
expected from l|2.10|l . corresponding to a higher string 
density. The rate of increase of the network length-scale 
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FIG. 3: The evolution of the typical string length-scale as a 
function of N. 



may now be extracted from the simulations by writing 
the late time behaviour of the simulation length-scale as 
£sim — oct, which can then compared to the asymptotic 
analytic form (|2.10() 



X 



1 - K N 



-f 



(3.11) 



In Fig. we plot a as a function of TV for the simula- 
tions given in Fig. |3 with the analytic form shown on 
the graph for comparison. Although the match is by 
no means perfect it does illustrate the general behaviour 
that £ increases more slowly as the number of string 
species increases, and also that the gradient a asymp- 
totes to a constant rather than zero, meaning that even 
for an arbitrarily large number of species the network 
will still scale rather than freeze. It is worth noting that 
our model underestimates the value of £tv, correspond- 
ing to an overestimate of the scaling density. This means 
that the network has another mechanism by which it can 
lose energy. One possibility, identified by j2]jj is that 
the network will lose energy from the binding energy of 
two string species as they merge into a bound state, this 
could account for our overestimate of the scaling density. 

We have repeated the simulations for a radiation domi- 
nated Universe and present the results in Fig. [5] Plotted 
here is the ratio of the energy density in the string net- 
work to that of the background radiation as a function 
of time. For each value of N there are two curves shown, 
each with a different initial value for the ratio of string 
density to radiation density. The figure shows that for a 
given N each curve seems to approach the same constant 
value, indicating that there is a unique scaling solution 
for each N. Unfortunately the simulations run out of 
dynamic range before we are able to get accurate values 
for the scaling density, however we can see from the plots 
that the scaling density increases with N. As the density 
in radiation scales as a -4 , with a ~ yt we have that the 
network density during scaling varies as p s ~ 1 /£ 2 ~ 1 /t 2 
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FIG. 4: The slope of f (t, N) as a function of N. 



which gives us the length-scale increasing as £ ~ at as 
before. Moreover, the simulations in the radiation back- 
ground indicate that the scaling density increases with N 
implying that a reduces as N increases, consistent with 
our analytic approach. 
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FIG. 5: Scaling in radiation as a function of N. 



We can provide further evidence for a scaling solution 
by measuring the total length of string contained in loops 
and the total length in links. In a scaling solution there 
should only be a single length-scale describing the prop- 
erties of the network, as such we would expect the ratio 
of the amount of string in loops to the amount of string 
in links to approach a constant. To test this we focused 
on the case of N = 3 and ran 100 simulations in order 
to minimize statistical errors. The results arc plotted in 
Fig. [S] where we see how the ratio of loop length to link 
length evolves in time, we can clearly see that the ratio 
tends to a constant, furnishing us with further evidence 
that the network approaches scaling. We now proceed to 
describe the relative length distributions of the loops and 
links in these non Abelian networks. 
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FIG. 6: The ratio of loop length to link length for N = 3. 



IV. LENGTH DISTRIBUTION 

We have presented evidence that the string networks 
enter a scaling regime which can be described by a single 
length-scale. This however does not mean that all the 
loops and links have the same length given by £, so we 
now turn our attention to the distribution of the lengths 
in loops and links. It is well known that the more usual 
Abelian cosmic string has a pow er law distribution of 
loop length, n(l)dl ~ l~ 5 / 2 dl [28j |. however our case for 
multip le sp ecies of string is more like the situation stud- 
ied in [3(j with three string species that could connect 
via heavy monopoles. There it was found that the distri- 
bution was more closely modelled by an exponential. Ex- 
ponential distributions for string lengths were also found 
in p9j where the statistics of a random distribution of Z3 
strings were studied. Here, we have focused on the case 
N = 3 and performed 100 simulations in order to get 
good statistics for the length distribution of loops and 
links. The results are given in Fig. Q where we have 
plotted lnn(i) vs I, when the system has approached the 
scaling regime. At large I the nearly-straight line of the 
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FIG. 7: The distribution of loop/link lengths for N = 3 



plot clearly shows that the distributions are more accu- 
rately described by exponentials rather than power laws. 
We may model this by an exponential, dimensionless loop 
production function f(l/£) by adapting the arguments 
found in 28] for the Abelian string. This function f(l/£) 
describes the loss from the network of strings of length 
I. We consider the network to be in the scaling regime 
described by the length-scale £ such that (|2.2|) holds for 
the network. The loop production function f(l/£) is de- 
fined such that nf(l/£)dl/l is the energy lost in loops of 
length between I and I + dl in volume £ 3 per unit time. 
This then means that the loop density evolves according 
to 



PL 



-3Hp L + gf if(l/0/Z 4 



(4.12) 



where g is a Lorentz factor to account for the loops cre- 
ated with non-zero centre of mass kinetic energy. We 
may rewrite this expression in terms of the dimension- 
less parameter x = l/£ and integrate to find 



pL(l,t) 



gill 



3f-3 



dx x 



2-3l/ 



/(*) 



(4.13) 



l /act 



Here we have taken a(t) ~ t v , £ = at, note that the 
radiation era corresponds to v = 4. In order to get an 
exponential distribution of loop lengths we must there- 
fore have a loop production function which contains an 
exponential, for example we may have 



f(x) — bx p exp(— fix), 



(4.14) 



where b, p and fi are constants. Taking this form and 
using n(l,t) — pz,(l,t)/l, l|4.13l) gives: 
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where in the final expression we have taken the limit 
x » 1. In principle we should be able to determine the 
dependence onp from the numerical simulations, but this 
is not possible as it stands due to the statistical errors 
in the data at large loop length. In other words we can 
not yet determine the form of any power law correction 
to the pure exponential distribution. 

We finally return to the distribution of the total string 
lengths in links and loops in the scaling regime, as shown 
in Fig. © . We can understand the late time behaviour of 
this Figure based on our analytic approximations. If we 
consider the link length, Li; n ks to be simply the difference 
between the total length of string L and the length of 
string in loops, L\ oops , then we can write : 
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(4.16) 



Now in a large volume V, where the string tension is \x 
we know that 
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where pi 00 ps is the energy density in loops defined from 
Eq. jOS) by 

pat 

Pioo ps = / p L (l,t)dl (4.18) 

Jto 

the limits of integration accounting for the largest loops 
excised at time t and the smallest loops excised at time 
to. Performing the integral we find 

bgfj, 1 

= aH 2 ^3-3,+/ (ft? i - f ) ( 4 ' 19 ) 

where 

F(fl,p,t)= f x 3{u - 1) r(3(i'-l)+p,l3x). (4.20) 

Now, F is independent of t at late times, hence we can 
extract the time dependence in Eq. I|4.19|l without know- 
ing the full form of F. Now given p = and f = at, 
we find the late time behaviour in Eq. 14.17( 1 

— — = =P. t (4.21) 

Sloops bgF 

which is independent of time, as required. 

V. CONCLUSIONS 

In this paper, in an attempt to model (p,q) strings aris- 
ing out of compactifications of type IIB string theory, 
we have studied how a non-Abelian network of strings 
evolves. Such networks contain multiple vertices where 
many different types of string join together. For example, 
Z3 strings have triple vertices where three different types 
of string join together. The result of these more compli- 
cated vertices is that non-Abelian networks have a richer 
variety of re-connections than their counterpart Abelian 
networks. Given this rich structure, from a cosmological 
point of view, we need to know whether such networks 
approach a scaling solution where they track the dom- 
inant background energy density or simply freeze, and 



stretch with the expansion of the universe leading to a 
string dominated Universe. By re-considering the field 
theory model of [2^ , which can be thought of as contain- 
ing N basic strings and ^N(N — 3) bound states to give 
a total of ^N(N — 1) different species of string, we have 
shown both numerically and analytically that we always 
approach a scaling solution for all N, in other words the 
system does not freeze. All that happens is that as N 
increases, so does the scaling density, asymptoting to a 
fixed value independent of N for large N. Moreover, we 
have seen firm evidence that the associated distribution 
of loops and links is not governed by a power law as 
is the case for Abelain and Z2 strings. Rather, as ini- 
tially observed on smaller lattices by Aryal et al |23| and 
Vachaspati and Vilenkin j^jj, the distribution drops off 
exponentially fast with loop size. We have attempted to 
explain this distribution based on our one-scale model, 
and argued why the loop-link ratio tends to a constant 
at late times for the case of Z3 strings. There is much 
that remains to be done. Perhaps the most obvious thing 
is to allow the strings to have different tensions in this 
field theory model. In a very nice paper, Tye et al have 
been doing just that in the context of modelling cosmic 
superstrings |27| . It would be worthwhile to see whether 
the field theory mimics their results or not. In particular 
once we have different tensions, it would be interesting 
to see how much string of each tension exists in a scaling 
network. One of the nice features about our approach is 
that we can produce the strings in an initial phase transi- 
ton and watch them evolve. It would be interesting to see 
how quickly the heavier tension strings appear, starting 
say from a network of light strings. 
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